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§ ■ ABSTRACT 
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o ! 

D , We consider a simple Newtonian model of a steady accretion disk around a 

Q ! black hole. The model is based on height-integrated hydrodynamic equations, 

! a-viscosity, and a pseudo-Newtonian potential which results in an innermost sta- 

I— i! ble circular orbit (ISCO) that closely approximates the one predicted by GR. We 

p ^! find that, as the disk thickness H/R or the value of a increases, the hydrody- 

Q ! namic model exhibits increasing deviations from the standard thin disk model 

^ \ of Shakura & Sunyaev. The latter is an analytical model in which the viscous 

\ torque is assumed to vanish at the ISCO. We consider the implications of the 

\ results for attempts to estimate black hole spin by using the standard disk model 

^ \ to fit continuum spectra of black hole accretion disks. We find that the error in 

\ the spin estimate is quite modest so long as H/R < 0.1 and a < 0.2. At worst 

^ I the error in the estimated value of the spin parameter is 0.1 for a non-spinning 

^ I black hole; the error is much less for a rapidly spinning hole. We also consider the 

^ I density and disk thickness contrast between the gas in the disk and that inside 

I the ISCO. The contrast needs to be large if black hole spin is to be successfully 

^ I estimated by fitting the relativistically-broadened X-ray line profile of fluorescent 

. !^ I iron emission from reflection off an accretion disk. In our hydrodynamic models, 

^ I the contrast in density and thickness is low when H/R > 0.1, suggesting that the 

I iron line technique may be most reliable in extremely thin disks. We caution that 

these results have been obtained with a viscous hydrodynamic model. While our 
results are likely to be qualitatively correct, quantitative estimates of, e.g., the 
magnitude of the error in the spin estimate, need to be confirmed with MHD 
simulations of radiatively cooled thin disks. 

Subject headings: X-ray: stars — binaries: close — accretion, accretion disks — 
black hole physics 
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1. INTRODUCTION 

Recently, we reported spin estimates of three black holes (BHs) in Galactic X-ray bi- 
naries (Shafee et al. 2006; McClintock et al. 2006; hereafter S06, M06). The results were 
obtained by fitting the soft X-ray continuum spectra of these systems in the thermal state 
(Remillard & McChntock 2006) to a general relativistic, multicolor blackbody, thin disk 
model {Kerrbb, Li et al. 2005), which includes the effect of spectral hardening (Davis et al. 
2005). In this method, which was pioneered by Zhang, Cui & Chen (1997), we assume a 
razor-thin disk that terminates at the innermost stable circular orbit (ISCO). In addition, we 
assume that the viscous torque vanishes at the ISCO and that there is no energy dissipation 
or angular momentum loss inside the ISCO. These are standard assumptions in the theory of 
accretion disks (e.g., Shakura & Sunyaev 1973; Frank, King & Raine 2002), and correspond 
to what we refer to in this paper as the "standard disk model." However, there has been 
debate in recent times as to the validity of the assumptions. 

The stress responsible for angular momentum transport in a thin accretion disk is likely 
to be magnetic (Balbus & Hawley 1991). If this is the case, an argument could be made 
for a non-zero stress at the ISCO, coupled with considerable dissipation near and inside the 
ISCO (Krolik 1999; Gammie 1999). These effects could cause important deviations from the 
standard disk model (Krolik & Hawley 2002), perhaps invalidating our spin determinations. 

Afshordi & Paczyhski (2003), following earlier work by Abramowicz & Kato (1989) 
and Paczyhski (2000), suggested that the torque at the ISCO increases with increasing disk 
thickness. Motivated by their work, we argued in M06 that deviations from the standard 
disk model are likely to be serious only for thick disks. We thus restricted our attention 
to relatively thin disks with height-to-radius ratios of H/R < 0.1. The present paper is an 
attempt to verify whether or not such thin disks do indeed behave like the standard disk 
model. 

In addition to the debate over the vahdity of using the standard disk theory to model the 
continuum spectra of realistic disks, another relevant issue in attempting to estimate BH spin 
is the relative merit of the continuum fitting method compared to fitting the relativistically- 
broadened fluorescent iron line in the X-ray spectrum. Both methods have been proposed as 
a means of estimating BH spins, and it is of interest to understand how well the assumptions 
of each are satisfied by real disks. The models currently used by the iron line method 
assume that the line emissivity peaks at the ISCO, drops abruptly to zero inside the ISCO, 
and decreases steeply as a broken power- law outside the ISCO (e.g., Brenneman & Reynolds 
2006, hereafter BR06). This requires, among other things, a significant drop in matter 
density (Fabian 2007) or disk thickness (Nayakshin et al. 2000, 2002) inside the ISCO. A 
second motivation for the present paper is therefore to check the validity of the assumed line 
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emissivity profile. 

Our analysis is based on a non-relativistic hydrodynamic model of an accretion disk. 
We present global numerical solutions of the differential equations governing the fluid flow, 
assuming that the accretion disk is steady, axisymmetric and in hydrostatic equilibrium in the 
vertical direction, and using a pseudo-Newtonian model for the gravitational potential. We 
do not include magnetic fields explicitly, but assume an effective viscosity described by the a 
prescription (Shakura & Sunyaev 1973). We also assume an adiabatic index 7 = 1.5, which 
corresponds to approximate equipartition between gas and magnetic pressure (Quataert & 
Narayan 2000). 

Our primary interest is in accretion disks in the rigorously defined thermal state (see 
Table 2 in Remillard & McClintock 2006) with HjR < 0.1, as these are the systems of 
most interest for our work on BH spin (M06). Since the value of the viscosity parameter 
a for such disks is a matter of debate, we try different constant values: a = 0.01, 0.1, 
0.2. We also consider a variable-a prescription (eq. 22) inspired by the MHD simulations 
of Hawley & Krolik (2002, hereafter HK02). For non-spinning BHs, we use the pseudo- 
Newtonian potential of Paczyhski & Wiita (1980; hereafter PW80), and for spinning black 
holes we use the pseudo-Kerr model of Mukhopadhyay (2002). The numerical framework for 
our calculations is similar to that used by Narayan, Kato & Honma (1997), viz., we use a 
relaxation method to solve the equations from the sonic radius Rs to the outer edge of the 
disk (~ 10^ Rs), and we then integrate inward from Rg to the event horizon. 

The paper is organized as follows. We discuss in §2 the theory and computational 
method. We then discuss in §3 our numerical disk solutions, focusing on the magnitude of 
the stress at the ISCO, the amount of viscous dissipation near and inside the ISCO, and the 
density and disk thickness contrast across the ISCO. We then compute in §4 the emitted 
spectra of our numerical disks for different values of H/R and a and investigate the error 
we make when we estimate the spin of a BH via the continuum fitting method assuming the 
standard disk model. We conclude in 55 with a discussion. 



2. THE MODEL 
2.1. Gravity 

In order to focus our attention on the key physics of the problem, and to avoid being 
distracted by technical details, we consider a simple viscous hydrodynamic accretion disk 
in a Newtonian gravitational potential. Since the presence of an ISCO is essential for our 
analysis, we simulate relativistic gravity in this Newtonian model by means of a modified 
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gravitational potential. For a non-spinning BH, we make use of the PW80 potential: 

GM 



R — 2i?g 



where M is the BH mass, G the gravitational constant and = GM/c^. The Keplerian 
angular velocity JIk at a radius R from the BH is 

In the case of a spinning BH we use the pseudo-Kerr model of Mukhopadhyay (2002) in 
which the gravitational acceleration of a test particle in a Keplerian orbit at a distance R 
from the BH is 

where r = R/Rg, a* = a/M = J/{GM'^/c) is the dimcnsionless spin of the BH, and —1 < 
a* < 1. The Keplerian angular velocity at radius R is then 



^ (r^-2a,^ + a;) 

GM r^{^{r - 2) + a.) " ^ ^ 



2.2. Hydrodynamics 

We assume a steady axisymmetric disk in hydrostatic equilibrium in the vertical di- 
rection. In the equations that follow, which have a long history in accretion disk theory 
(e.g., Paczyhski & Bisnovatyi-Kogan 1981; Muchotrzeb & Paczyhski 1982; Kato, Honma & 
Matsumoto 1988; Abramowicz et al. 1988; Popham & Narayan 1991; Narayan & Popham 
1993; Chen & Taam 1993; Narayan et al. 1997; Chen, Abramowicz & Lasota 1997), we 
denote density, sound speed, radial velocity, angular velocity, Keplerian angular velocity, 
and vertical half-thickness by p, Cg, vr, fl, Qk, and H, respectively. All these parameters are 
taken to be functions of the cylindrical radius R only. Because of the assumption of steady 
state, the Lagrangian time derivative D/ Dt — d/dt + v -V becomes D/Dt — v-^d/dR. After 
vertical and then radial integration the continuity equation takes the form: 

Anpv-s^RH = —M = constant , (5) 

where H — Cs/JIk- The momentum equation is 



p(v • V)v = -VP - pV$ + pn^R + pV -a , 



(6) 
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where a is the stress tensor. We assume that the only non-zero component of a is (Tr$ = —aP 
(a prescription, Shakura & Sunyaev 1973), where P is the total pressure and we write 
P — pCs2. The radial component of the momentum equation gives 

and conservation of angular momentum gives 

R dR^ ' Rm dR ^ ' 

The latter equation can be integrated to obtain 

nR'-j = -'^, (9) 

where VLR^ is the specific angular momentum of the gas at radius R and j is an integration 
constant. We can interpret j as the specific angular momentum of the accreting gas at the 
radius where the stress goes to zero. 

Lastly, we write the energy conservation equation in terms of the Lagrangian derivative 
of the specific entropy, 

pT^ = q+-q-= fq+. (10) 

Here s is the specific entropy per unit mass, and g"*" and q~ are the volume rate of heating 
and cooling of the gas, respectively. Following Narayan et al. (1997) wc take the cooling 
rate to be a factor (1 — /) of the heating rate. Narayan et al. used / = 1 because they 
were modeling advection-dominated accretion flows. Since we are interested primarily in 
thin disks, we use small values of /, i.e., substantial coohng, and we tune the value of / 
to achieve the desired disk thickness (eq. 19). The heating of the gas is due to viscous 
dissipation, which gives q'^ — vaRdil/dR. Using the relationship e — P/{'j — 1), where e is 
the thermal energy per unit volume and 7 is the adiabatic index (we use 7 = 1.5), we can 
write 

Ds pvR d4 2 dp 

pi — = c:vo — . (11) 

Dt 7-1 dt " dR ^ ' 

Thus, the energy equation takes the form 

pvK del 2 dp ^ dVt 

- '^sVr-tb = -/Wi?— . (12) 



j-ldR ' ^dR ' ^ ' dR 
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2.3. Boundary Conditions and Numerical Method 

We use a relaxation method to obtain numerical solutions of the above differential 
equations. In the computations, we define x = R/Rs as the spatial variable and covered 
the region ,t = 1 to 10^ using 1000 grid points. The grid has a non-uniform spacing, 
with more grid points near the inner boundary x — 1. In solving the equations, we set 
—ATTpv^RH — M — 1, and in order to simplify the equations we substitute for p using 
equation (5). Thus we are left with three unknown functions of R: i'r(-R), cl{R), and 
Q{R). In addition, we have two unknown constants, j and Rg, which we treat as eigenvalues. 
To solve for these quantities, we use equations (7), (9), and (12), supplemented with five 
boundary conditions. 

Narayan et al. (1997) showed that solutions of the disk model described in §2.2 tend 
to be nearly self-similar over a wide range of radius. Assuming self-similarity (following 
Narayan & Yi 1994), we can obtain the following analytic solution of the equations (the 
subscript "SS" refers to self-similar): 



^ iR)-^™ 2 ^- 5/3-7 




vk,ss{R) = VQ\j , vq = -a\/ ^ , (14) 



We use this self-similar solution to set boundary conditions at the outer boundary i^out^ 



^^R(^out) = VqJ^^ , (16) 
V -Kout 

Cs^(^out) = , (17) 

n{Rout) = ^^O^^K . (18) 

Prom the above relations it can be shown that, at the outer boundary, the vertical scale- 
height H satisfies 

H 



R Y 5 + 2e' + aVe' ■ ^^^^ 

Therefore, for a given value of 7, we can vary the disk thickness H/R by changing /. For 
7 = 1.5, /= 0.000035 and 0.0035 give H/R = 0.01 and 0.1, respectively Once set at the 
outer edge, the value of H/ R remains constant over most of the disk, becoming smaller only 
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near and inside the ISCO. Note that, for the thin disk models that we consider in this paper 
which have H/R < 0.1, the advection parameter / is very much less than unity. This means 
that radiative cooling (which is oc 1 — /) dominates by a huge factor over energy advection 
(oc /). We briefly discuss thicker advection-dominated solutions in §5. 

The inner boundary is at the sonic radius, R — Rs, which is a singular point of the differ- 
ential equations. Following standard methods, we obtain the following regularity conditions 
at Rg'. 

"l - = . (20) 

^ ^ ^ '7 + 1 dR J '7 + 1 di? ^ ^ 

Equations (16)-(18), (20)-(21) provide the five boundary conditions we need to find a unique 
solution. Once we have obtained the solution between R — Rs and R — Rout via the 
relaxation method, we use the solution at R — Rg ets initial conditions and integrate the 
equations from Rg down close to the BH event horizon. 

We should emphasize that we do not set any boundary condition at the ISCO. Instead, 
we apply the boundary conditions at the sonic radius, whose position is computed self- 
consistently for each solution. Further, even at the sonic radius, the viscous torque is not 
set to zero — the torque is computed self-consistently and is allowed to continue smoothly 
inside the ISCO. The numerical solutions wc obtain arc thus superior to the standard disk 
model and can be used to check the validity of the latter. In particular, we can estimate 
what error one makes in the standard disk model as a result of the zero-torque boundary 
condition. 



3. RESULTS 

3.1. Numerical Solutions 

Figure 1 shows model results for a non-spinning BH. We consider two disk thicknesses: 
H/R — 0.01 (sohd lines), and H/R — 0.1 (dotted lines). In all four panels the vertical fine 
shows the position of the ISCO {R — 6i?g). We use G = M = c = 1, so that the unit 
of velocity and time are c and GM/c^, respectively, and set M = 1. Most of our models 
correspond to a constant value of a. However, we also consider a model in which a varies 
as a function of R, 

16.8 , , 
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which closely reproduces the effective profile of a found by HK02 (see their Fig. 4) . We refer 
to this as the "variable-a model." 

Figure la shows the variation of the sound speed squared as a function of radius R. For 
a given thickness, the different a models overlap at large radii and are only distinguishable 
in the inner region of the disk. Here and in the figures that follow, the magenta, blue, red 
and green lines refer to the a = 0.01, 0.1, 0.2 and variable-a models, respectively. Figure lb 
shows the radial infall velocity of the accreting gas. We see that, between the ISCO and 
the event horizon, vr increases rapidly regardless of the value of a. The variable-^ model 
almost completely overlaps with the a= 0.1 model even at large radii. Figure Ic shows the 
angular velocity Q and Keplerian angular velocity Qk- The profiles of Q for the different 
values of a and H/ R are not distinct and are represented by the single dotted line. The sohd 
red line corresponds to the Keplerian angular velocity. Note that the gas orbits in a nearly 
Keplerian fashion until it reaches the ISCO. Thereafter, the hydrodynamic forces maintain 
an orbital motion that becomes increasingly sub-Keplerian as the gas approaches the event 
horizon. Figure Id shows the gas density p as a function of radius. As in the case of the 
sound speed (Fig. la), the density reaches a maximum outside the ISCO and then decreases 
rapidly near the event horizon. In this plot, too, the variable-o; model coincides with the 
Q! = 0.1 model. 

Figure 2 is in the same format as Figure 1 and presents our results for a spinning BH 
with = 0.95. The principal difference from the previous figure is that the ISCO (vertical 
dashed line) is now located at i? = 1.937i?g. We consider the same values of a and H/R as 
in Figure 1, but there is no variable-a model in this case because HK02 considered only a 
non-spinning BH. 

3.2. Matter Density, Disk Thickness and the Iron Line Method 

Before presenting our main results in the following subsections, we briefiy consider the 
implications of our models for the determination of spin via the iron line method. The 
source geometry and illumination law for producing the fluorescence iron line are probably 
the largest uncertainties in the hue fitting method (Reynolds & Begelman 1997). If we assume 
the steepest law that is suggested by Reynolds & Begelman (1997), then the irradiating flux 
Fx ~ R"^. Let us write the emissivity function in the form /pe-^x, where /pe is an efficiency 
factor. In this section we investigate if the existing models of /pe in the literature agree with 
our hydrostatic models. 

The currently favored iron fine models (BR06) assume that the iron line emission is 
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restricted between i?isco and an outer radius i?out and that, within this region, the hne 
profile is fitted by a broken power law. BR06 find that the emissivity varies as ~ between 
the break radius i?br and i?iscoi a^nd as ~ between Ri,^ and i^out- For Fx ~ R~^, this 
implies the following form for the efficiency function: 

/Fe(BR06) = , if < -RiSCO , 

= 1/R^ , if i^isco <R<S Risco , (23) 
= constant , if 3 Risco < R ■ 

Below we discuss the two main theories regarding the physical parameters that might affect 
the emissivity profile. 

Constant density models (Ross & Fabian 2003; Zycki et al. 1994; Ross, Fabian & Young 
1999) predict that the line emissivity is dependent on the ionization parameter, which is 
proportional to Fx/ p, where p is the gas density and Fx is the illuminating flux. It is 
argued that the gas density drops to very low values inside the ISCO. As a result, the region 
inside the ISCO has a very high ionization parameter, which in turn produces negligible 
iron line emission (Reynolds & Begelman 1997; Young et al. 1998; Fabian 2006). In this 
case, one would expect /pe to be inversely related to ionization, i.e., /pc should be a function 
of p{R)/ Fx{R) oc p{R)R^ . More detailed calculations that solve for the vertical structure 
of the disk under hydrostatic equilibrium (e.g. Nayakshin et al. 2000, 2002) suggest that 
the line emission depends on a "gravity factor" ~ [H / R^)Fx- If that is the case, then for 
Fx ~ R~^-i one expects the efficiency function /pe to be proportional to H. 

In Figure 3, we compare the BR06 efficiency function /fc(bro6) (e<l- 23) to those sug- 
gested by our hydrostatic models, in the context of the constant density and gravity theories 
mentioned above. Figure 3a shows p{R)R^ as a function of radius. The hne types/colors for 
the various models are the same as those defined in Figure 1. Superimposed on our density 
profiles is a thick short-dashed black hne that represents /Fe(BR06)- For both H/R — 0.01 
(solid hues) and H/R — 0.1 (dotted hues), and all values of a, we note that p{R)R^ is an 
increasing function of radius, implying that /pe should also increase with increasing radius. 
There is no apparent reason why /pe should increase so steeply near the ISCO, or decrease 
at large radii, as suggested by /pe(BR06)- 

Figure 36 shows a similar plot for a rapidly spinning BH with a* = 0.95. We notice the 
same trends as in Figure 3a. In this case, we also notice that (especially for H/R — 0.1), 
instead of becoming negligible at the ISCO, p{R)R? decreases gradually as one passes the 
ISCO and moves closer to the event horizon. Therefore, one does not expect /pe to drop 
abruptly to zero at the ISCO. 

In Figures 3c and 3d, we consider the disk thickness in the inner region. In our models. 
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the disk has a more or less constant thickness specified by H/R outside ~ 100i?g, and we 
vary this "outer thickness" by changing the value of / (eq. 19). However, in the inner region, 
the disk gets thinner. Figure 3c shows H as a, function of R for a non-spinning BH. The 
top panel shows a disk with outer thickness of 0.01 and the bottom one shows a disk with 
outer thickness of 0.1 for the choices of a specified in §3.1. In the thinner case, there is 
an abrupt drop in H, which would likely quench the iron emission from inside the ISCO. 
For the thicker case, however, the value of H decreases gradually and remains significant far 
inside the ISCO at 3it!g. Thus, these models indicate that the region within the ISCO may 
contribute a significant fraction of the total iron line emission and, also that it is difficult to 
justify the steeply falling form of /Fe{BR06) ■ As shown in Figure 3d, the results for a BH with 
a* = 0.95 (-Risco = 1.973i?g) are very similar. Again, for H/R — 0.1 the disk thickness H 
decreases gradually near and within the ISCO. 

We hasten to add that this is a very simple model of an accretion disk, perhaps too 
simple to address "surface phenomena" such as fiuorescent iron line emission. Modulo this 
important caveat it seems that, for reasonable values of the model parameters, the iron line 
emission does not necessarily end at the ISCO, nor does it vary with radius outside the ISCO 
with anything like the functional form assumed in current fits of iron line data (e.g., BR06). 

3.3. Viscous Stress Near the ISCO 

Figure 4 shows the vertically integrated stress 2HaP for a non-spinning BH. As shown 
in Figure 4a, all the models corresponding to a very thin disk are in close agreement with the 
standard model, i.e., the stress nearly vanishes at the ISCO even though we do not require 
this of the model. For the thicker disk shown in Figure 46, the stress near and inside the 
ISCO increases, the effect becoming more important for larger values of a. Interestingly, for 
a — 0.01, the magnitude of the peak stress is actually smaller than that predicted by the 
standard disk model. 

As shown in Figure 5, our models for a spinning black hole display essentially this same 
dependence of stress on H/R and a. In both Figures 4 and 5, the presence of a non-zero 
viscous stress inside the ISCO imphes a contribution to the observed spectrum that is not 
accounted for in the standard disk model. In §4 we investigate the magnitude of this effect. 

We now consider the effect of a and disk thickness on the eigenvalue j (§2.2), which is 
the specific angular momentum delivered to the black hole by the infalling matter. In the 
standard disk model, j is the Keplerian specific angular momentum at the ISCO because 
(i) matter is assumed to orbit at the Keplerian velocity and (ii) the stress is assumed to 
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vanish inside the ISCO. Neither assumption is made in our hydrodynamic models, and it is 
therefore of interest to consider how much the calculated values of j differ from the standard 
value. Table 1 summarizes the values of j for our different models. We note that j decreases 
with increasing H/R and a. That is, as the disk gets thicker or as a increases, more angular 
momentum is removed before matter falls into the BH. However, the effects are quite small, 
and the deviations are less than 1% in all cases. 



3.4. Dissipation Inside The ISCO 

In the previous section we showed that the stress at the ISCO is small, but non-zero, 
and that it increases with disk thickness and a. We now consider the energy dissipation 
profiles of our model disks for different values of a and H/ R. Figure 6 shows the quantity, 

dR d\n{R) dR ^ ^ ' 

as a function of R. Here L is the luminosity and D{R) the energy dissipated per unit 
time per unit surface area of the disk. Figures 6a and 66 show RdL/dR vs it! for a* = 0, 
while Figures 6c and 6d show the results for a* = 0.95. The sohd black lines show the 
standard disk model with zero torque at the ISCO. For the thin disk with H/R — 0.01 and 
for all values of a, our models are indistinguishable from the standard model, which thus 
provides an excellent description of the flow in this case. However, for the thicker disk with 
H/R = 0.1, our numerical models deviate somewhat from the standard disk model. We note 
in particular that larger values of a arc associated with more dissipation near the ISCO and 
larger deviations from the standard disk model. 

In Table 1 we summarize the total luminosities of the different models for a given mass 
accretion rate M . We note that none of the luminosities of our models deviates by more 
than 4% from that of the standard model. 
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Luminosities and Angular Momentum Eigenvalues of the Numerical Disk Models 



a, H/R 


a 


-^TOTAL 

(Mc^) 


-f'TOTAL(STD) J 

(Mc^) 


isTD 


0.01 


0.01 


0.0624 


0.0625 3.6744 


3.6742 




0.1 


0.0625 


3.6735 






0.2 


0.0626 


3.6727 






\'arial)le 


0.0626 


3.6730 




0.1 


0.01 


0.0610 


3.6839 






0.1 


0.0633 


3.6609 






0.2 


0.0650 


3.6456 






variable 


0.0646 


3.6518 




0.95 0.1 


0.01 


0.2091 


0.2144 2.3372 


2.3311 




0.1 


0.2171 


2.3237 






0.2 


0.2227 


2.3146 





The subscript STD refers to the standard thin disk model. 



4. DISK SPECTRA AND THE EFFECT ON BH SPIN ESTIMATION 

In the standard disk model, the viscous dissipation is assumed to vanish at the ISCO. 
As a result, the emitted flux also vanishes at the ISCO, and no radiation is emitted from the 
region of the flow between the ISCO and the event horizon. For a given BH mass, the radius 
of the ISCO is a well-known and monotonically decreasing function of a*, e.g., for a^— 0, 1, 
the ISCO is located at 6i?g, li?g, respectively. As discussed in Zhang et al. (1997), S06 and 
M06, the radius it!in of the inner edge of the disk can be estimated from observations. For a 
BH of known mass, this radius can be expressed in units of R^. and if the disk inner edge is 
located at the ISCO, then R^^jRg determines the spin parameter a*. 

From the calculations presented in this paper, we sec that for a very thin disk [Hj R = 
0.01) the viscous dissipation does indeed become negligible inside the ISCO and the dissipa- 
tion profile RdL/dR is identical to that predicted by the standard disk model. Thus for such 
systems we expect our estimates of BH spin to be quite accurate. However, we do notice a 
difference for thicker disks with say H/R ~ 0.1. Using the standard disk model to fit the 
observed spectra of these systems will lead to an error in our estimate of the radius of the 
ISCO. We now try to quantify this error. 
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For each of our disk solutions, we have calculated the emitted spectrum assuming that 
the disk emits like a blackbody at each radius. The temperature profile T{R) of the disk 
surface can be calculated from dL/dR using: 

{l-f)^ = ^7raRT\R), (25) 

where a is the Stefan-Boltzmann constant. This can be used to calculate the observed 
spectrum of the disk by integrating over the entire disk: 

_27rcosi 2hu^RdR 

- ]J2 ^2^{hv/kT{R))-l) ' V26j 

^ -dinner 

where h is the Planck constant, c the speed of light, k the Boltzman constant, D the distance, 
i the angle of inclination, v the frequency, and i^mner the radius of the inner boundary of the 
disk, near the event horizon. 

Figm^c 7 shows our calculated spectra for a BH with mass M — lOM© and distance 
D = 10 kpc. In each panel the solid curve shows the spectrum from a standard disk model 
with the appropriate pseudo-Newtonian potential. As before, we have considered three 
constant values of a: 0.01, 0.1 and 0.2, for both the spinning and non-spinning cases, and 
an additional variable-a model for the non-spinning case. 

Figures 7a and 76 show the calculated spectra for the case of a non-spinning BH. For 
H/R = 0.01, we see that the calculated spectra for all four models of a overlap with the 
spectrum calculated via the standard disk model. Therefore, we can conchidc that, for such 
very thin disks, the standard disk model is a very good approximation and that the choice 
of a cannot be a major source of error in estimating BH spin. The different curves are more 
distinct in the case of a thicker disk with H/R — 0.1 (magenta, blue, red and green hues 
show the a — 0.01, 0.1, 0.2, and variable-a models). The differences are especially noticeable 
at high photon energies, where larger values of a give higher fluxes. Figures 7c and 7d show 
H/R — 0.1 disk spectra for a* = 0.8 and 0.95. 

To estimate how the spectral distortions might affect BH spin determination, we pro- 
duced spectral data files for our models using an RXTE response file and analyzed the data 
with XSPEC version 12.2.0. These "fake" data files were fitted with the XSPEC model 
Diskpn (Gierlinski et al. 1999), which uses the standard disk model with the PW80 poten- 
tial and a zero torque boundary condition at the ISCO. Diskpn has three fit parameters: 
Tmax, Rin/Rg and normalization K = M'^ cos i / D'^ (3 , where M is the mass, D the distance, 
i the angle of inclination, and f3 the color correction factor. We are interested in the case 
when the inner edge of the disk coincides with the ISCO. Thus, since Diskpn considers a 
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non-spinning BH, we set Rm — 6Rg. The constant K can then be rewritten as 

V8.86 X 106cmy VlOkpcy (3 ' ^ ' 

In the above expression, 8.86 x 10^ cm corresponds to 6i?g — i?isco for a non-spinning black 
hole with M = IOMq. From the value of K obtained from spectral fitting, one can calculate 
for each model using equation (27). Using this value of one can then calculate the 
BH spin for which the ISCO would be located at that radius. This is the spin that one 
infers from the fake spectral data, under the assumption that the standard disk model is 
correct. Since the model was calculated with full viscous hydrodynamics as described in 
earlier sections, the spin value derived assuming the standard disk model will be different 
from the true BH spin (a* = in this case) . The difference between the two values represents 
the error in the spin estimate, Aa*, caused by our use of the simplified standard disk model. 

The results of this analysis are summarized in Table 2. The results correspond to 
M = IOMq, 10 kpc, cosi = 1 and /3 = 1. 



Table 2 



Errors in Spin Estimation of a Non-spinning BH 



H/R 


a 


i?in/(8.86 X 10^ cm) 


Aa* 


0.01 


0.01 


1.011 


-0.011 




0.1 


1.008 


-0.008 




0.2 


1.006 


-0.006 




variable 


1.007 


-0.004 


0.1 


0.01 


1.021 


-0.019 




0.1 


0.960 


0.037 




0.2 


0.920 


0.074 




variable 


0.879 


0.060 



Figure 8a shows the results in more detail for H/R = 0.01, 0.02, 0.04, 0.06, 0.08 and 
0.1. Figure 86 shows the results for the variable-a model as a function of disk thickness. We 
see that the error is larger for thicker disks and also for larger values of a. However, even for 
the thickest case we considered, H/R = 0.1, and the largest value of a = 0.2, the BH spin 
is overestimated by less than 0.1. Thus, in the case of a non-spinning BH the error is quite 
modest when one considers, for example, that both the radius of the ISCO and the binding 
energy at the ISCO differ only slightly (by 6%) for a BH with a* = 0.1. 
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Though a similar standard pseudo-Kerr XSPEC model is not available for fitting our 
spinning BH model spectra, it is still possible to estimate Aa^, by calculating the model 
luminosities. Gierlinski et al. (1999) showed that, for a non-spinning BH with the PW80 
potential, one can write: 

L = ^Mc^ = 35.7^i?fXa. , (28) 

where L is the luminosity, cr the Stefan-Boltzmann constant, (5 the color correction factor, R^^ 
the radius of the inner edge of the disk, and T^ax the peak temperature of the disk. Therefore, 
instead of calculating the multicolor blackbody spectrum of our models and fitting them 
with XSPEC, we could simply compare each hydrodynamic model with the corresponding 
standard disk model with the same Tmax- This gives the following estimate for the effective 
disk inner radius i?in of any given hydrodynamic model, 

-n-in - 7 -n^lSCO > i^yj 

-^standard disk 

where Lmodei is the luminosity of the model, i^gtandard disk is the luminosity of the standard 
disk with the same value of Tmax, and -Risco is the radius of the ISCO. The value of Ri^ 
obtained using equation (29) may then be used to calculate Aa*, as before. 

Figure 9a shows Aa* values calculated using both the full spectral fitting method via 
equation (27) and the simpler luminosity-temperature method described by equation (29). 
We see that the results are very close, indicating that the second method is a good proxy 
for the more detailed spectral method. 

For a spinning BH, equation (28) can be generalized to 

71(7 

L = eMc' ^ Co^I^A^ , (30) 

where e is the spin-dependent efficiency of the BH, and Cq is a constant. Therefore, equation 
(29) can again be used to estimate the effective R^^ and this can be used to obtain an estimate 
of the BH spin. 

Figure 96 shows Aa* for spinning BHs using this method. We show results for a* = 0.7, 
0.8, 0.9 and 0.95, and H/R — 0.1. For a given disk thickness and a, we see that the error 
in the spin estimate becomes smaller as the spin of the BH increases. For a* = 0.95, the 
maximum error is only ~ 0.01. 

5. DISCUSSION 

In this paper we studied the properties of a simple hydrodynamic model of an accretion 
disk using the a prescription for viscosity. We considered models with finite thicknesses H/ R 
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and different values of a. Our aim was to investigate how much the hydrodynamic models 
of thin disks deviate from the idealized "standard disk model" which assumes a vanishing 
torque at the innermost stable circular orbit (ISCO). 

We find that the deviations of the viscous hydrodynamic models from the standard disk 
model increase with increasing H/R and increasing a. However, even for H/R — 0.1 and 
a = 0.2, the largest values we tried for our thin disk calculations, the deviations remain 
modest. This is illustrated in Figures 4 and 5, which show how the stress profile deviates 
from that of the idealized standard disk model, and also in Figure 6, which compares the 
profiles of the viscous energy dissipation rate RdL/dR, Figure 7, which shows the multicolor 
blackbody spectra of the models, and Table 1, which gives some quantitative results. In all 
cases, we see that the detailed hydrodynamic models match the standard disk model quite 
closely. 

We were motivated to do this study because we and others have used the standard disk 
model to fit the continuum spectra of BH X-ray binaries in the thermal state in order to 
estimate the spins of the BHs. How much error do we expect in the estimated spin values as 
a result of the fact that a real disk deviates from the standard disk model? At least for the 
simple hydrodynamic models we have considered in this paper, the answer is that the errors 
are quite modest. 

Quantitative results are given in Table 2 and Figures 8 and 9. The error Aa* in the 
derived estimate of BH spin is at most ~ 0.1 in the case of a non-spinning BH and is much 
less for rapidly spinning BHs. These errors are comparable to or smaller than the errors that 
arise from uncertainties in our estimates of mass, distance and disk inclination (S06, M06). 

While these results are very encouraging for our program to estimate BH spin through 
fitting the continuum spectra of BH accretion disks in the thermal state, we must note 
some caveats. First and foremost, we have considered a highly simplified toy hydrodynamic 
model with an a prescription for viscosity. Real disks doubtless have magnetic fields, and 
the stresses associated with these fields probably do not behave like microscopic viscosity. 
Indeed, it is precisely this argument that has been used by Krolik (1999), Gammie (1999) 
and HK02 to question the zero-torque boundary condition at the ISCO. On the other hand, 
Paczynski (2000) makes an equally persuasive argument (based on the angular momentum 
conservation equation) that, so long as the shear stress is smaller than the pressure, a thin 
disk will always satisfy the zero-torque condition. 

In an attempt to include some of the effects of magnetic fields, we have considered a 
model in which we allowed a to vary with radius (sec cq. 22) in such a manner as to closely 
mimic the effective a obtained by HK02 from their MHD simulations. Even though in this 
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model a increases rapidly with decreasing radius, especially inside the ISCO, we found that 
none of our results changed. Based on this finding we cautiously suggest that the inclusion 
of magnetic fields may not significantly alter our conclusions. 

One question that needs to be addressed is why our results differ so much from those 
obtained by HK02. Prom MHD simulations of magnetized gas accreting in a PW80 potential, 
those authors concluded that the vertically integrated magnetic stress increases monotoni- 
cally with decreasing radius all the way through and inside the ISCO. This is dramatically 
different from the behavior we find, as a comparison of HK02's Fig. 10 with our Fig. 4 shows. 
A likely explanation is that we have limited our study to ihin disks {H / R = 0.01, 0.1) in 
which we simulated strong cooling by choosing a small value for the advection parameter / 
(see the discussion below eq. 19). HK02, by contrast, had no cooling in their MHD simula- 
tion, so their gas retained whatever energy was generated through shocks, making their disk 
thicker. 

In order to verify that this difference is important, we calculated models with larger 
values of / using our viscous hydrodynamic code. It is hard to know what effective value 
of / is most appropriate to match the HK02 simulation. Nominally, theirs was a fully 
advection-dominated accretion flow, since they had no cooling at all; this means that their 
simulation corresponded to / = 1. However, we do not know how well their code conserved 
energy. Therefore, we calculated three models with / = 1, 0.5 and 0.25, all with the variable 
a prescription (eq. 22) which most closely matches their stress profile. Figure 10 shows 
the resulting stress profiles. We sec that these advection-dominated models do exhibit a 
monotonically increasing stress inward, exactly as found by HK02 (their Fig. 10). The 
stress profiles are very different from those we find for cooling-dominated thin disks (our 
Figs. 4 and 5). Thus, we tentatively suggest that a large part of the difference between the 
results we find in this paper and those obtained by HK02 is related to the differing treatments 
of the energy equation of the gas, viz., cooling-dominated thin disk regime versus advection- 
dominated thick disk regime. In other words, we confirm the original insight of Abramowicz 
& Kato (1989), Paczynski (2000) and Afshordi & Paczynski (2003) on the strong relation 
between disk thickness and the stress at the ISCO. However, only a detailed MHD study of 
an accretion disk with significant cooling can tell for sure if this interpretation is correct, 
and to our knowledge nobody has carried out such a study. 

Another limitation in our work is that we used a Newtonian model and we simplified 
the thermodynamics of the gas in the disk via the advection parameter / (see eq. 10). 
However, doing the calculations in general relativity with full radiation thermodynamics 
will, we believe, introduce modifications only of order unity. The changes will be larger 
for a spinning BH, which we modeled with the Mukhopadhyay (2002) model, compared to 
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a non-spinning hole (PW80 potential), but we think the error will still be only of order 
unity. Therefore, calculating these effects in more detail will not greatly alter our qualitative 
conclusion that the standard disk model is adequate so long as the disk geometrically thin. 
Nevertheless, it would be useful to extend this work using a more complete set of disk 
equations, such as those employed in the study of slim disks (Abramowicz et al. 1988), and 
with the inclusion of general relativity (e.g., Abramowicz, Lanza & Pcrcival 1997). 

Note that we employed the two pseudo-Newtonian potentials mentioned above in the 
work reported here merely to obtain a ballpark estimate of the error associated with the 
zero-torque approximation. When we actually fit data to estimate the spin parameters of 
BHs (e.g., S06, M06), we use a detailed model (Li et al. 2005) which assumes the Kerr 
metric and includes all special relativistic and general relativistic effects. 

In our work on BH spin (S06, M06), we limited ourselves to disks with luminosities less 
than 30% of Eddington, which corresponds to vertical thicknesses H/ R < 0.1. The present 
study shows that this was a reasonable choice. For H/R < 0.1, the effects of gas physics 
and finite vertical thickness in our hydrodynamic models are not serious. Equally clearly, 
for thicker disks with H/R much greater than 0.1, the effects will be large; e.g., see Figure 
10. Therefore, one should be cautious about applying the standard disk model to disks 
more luminous than 30% of Eddington. For this reason, we believe the results obtained by 
Middleton et al. (2006) for the spin of the microquasar GRS 1915-1-105 should be taken with 
caution. 

Strong observational evidence that fitting the X-ray continuum is a promising way to 
estimate black hole spin comes from a long history of fitting the broadband spectra of black 
hole transients using the simple non-relativistic multicolor disk model (Mitsuda et al. 1984; 
Makishima et al. 1986), which returns the temperature Tin at the inner-disk radius In 
their classic review, Tanaka & Lewin (1995) give examples of the steady decay (by factors 
of 10-100) of the thermal flux of transient sources during which R^ remains constant. They 
remark that the constancy of R^ suggests that it is related to the radius of the ISCO. More 
recently, this evidence for a constant inner radius in the thermal state has been presented 
for a number of sources via plots showing that the bolometric luminosity of the thermal 
component is approximately proportional to (Kubota & Makishima 2001; Kubota & 
Makishima 2004; Abe et al. 2005; McChntock et al. 2007). In short, these non-relativistic 
analyses, which ignore spectral hardening (Davis et al. 2006), provide evidence for the 
presence of a stable radius, although they obviously cannot provide a secure value for the 
radius of the ISCO or even establish that the stable radius is the ISCO. 

We now consider the iron-line method of estimating spin. In this method, it is assumed 
that the line emission ceases abruptly at the ISCO, so an important question is whether or 
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not the gas inside the ISCO will fluoresce (Reynolds & Begelman 1997). The possibility of 
line emission from inside the ISCO is usually discounted on the grounds that the density will 
fall suddenly inside the ISCO, thus causing a sudden increase in the ionization parameter 
(Fabian 2007, and references therein). Alternatively, and to the same effect, it is argued 
that emissivity is related to the "gravity parameter" (Nayakshin 2000, 2002), and should 
depend on H. We see in Figure 3a that the density dependent function p{R)R^ does become 
neghgible inside the ISCO for the non-spinning BH. However, that is not the case for a fast 
spinning BH with a* = 0.95 (Figure 36) even for a small disk thickness of H/R — 0.1. Also, 
the radial dependence of H shown in Figures 3c and 3d implies that there should be emission 
from the inner region unless the disk is very thin {H/R — 0.01). 

An additional complication for iron line modeling is that the emissivity is assumed to 
vary as a broken power- law, with the maximum emission occurring exactly at the ISCO (e.g., 
BR06). Looking at Figure 3, such an ad hoc model would be hard to justify if the emissivity 
has anything to do with gas density or disk thickness. In contrast, the continuum-fitting 
model has the merit that it makes use of a physically motivated profile of disk emission 
RdL/dR which can be calculated from first principles in the standard disk model and which 
continues to be valid even in the more general hydrodynamic models described in this paper 
(Figs. 6, 7). 

This paper has focused on only one aspect of BH spin estimation, viz., the validity of 
assumptions made in various methods of spin determination regarding the hydrodynamical 
properties of the accretion disk. Of course, a successful determination of spin needs more 
than a valid disk model. It also requires high quality data and accurate determination of 
secondary system parameters. A discussion of these issues is beyond the scope of this paper, 
and the reader is referred to appropriate papers in the literature (e.g., M06; BR06). 

The authors thank Niayesh Afshordi and Jonathan McKinney for discussions and use- 
ful suggestions. We dedicate this paper to Bohdan Paczynski for his amazing insights in 
accretion theory and in numerous other areas of astrophysics. 
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Fig. 1. — Disk parameters for a non-spinning BH: (a) sound speed, (6) radial velocity, (c) 
angular velocity, (d) density. In panels a, b, and d, the solid and dotted lines correspond to 
H/R = 0.01 and H/R = 0.1, respectively; where distinct, the magenta, blue, red and green 
lines represent a = 0.01 , 0.1 and 0.2, and variable a (eq. 22), respectively. In all three 
panels, the variable-a model is nearly coincident with the a = 0.1 model. In panel c, the 
Keplerian velocity is plotted as a solid red line. Because the angular-velocity profile of the 
four models are nearly identical, we represent them by a single dotted line. The radius of 
the ISCO, R = QRg, is indicated in all four panels by the vertical dashed line. All numerical 
values correspond toG' = c = M = l,M = l. In panel c, the unit of angular velocity is 
{GM/c^)-\ In panel d, the unit of p is cy{G^M^). 
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Fig. 2. — Similar to Figure 1, but for a spinning BH with a* = 0.95. The vertical dashed 
line shows the ISCO at R= 1.937-Rg. There is no variable-a model in this case (see text), 
and hence the green line is not present in these plots. 
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Fig. 3. — Profiles of p{R)B? (panels a and h) and disk thickness H (panels c and d) in the 
inner regions. For the line types/colors defining the various models and the location of the 
ISCO, see Figures 1 and 2. Superimposed on our models in all of the panels is a thick black 
line that schematically represents the emissivity profile efficiency, /pe, assumed in the iron 
line work in which the emissivity cuts off abruptly inside the ISCO and falls off as a steep 
power law outside the ISCO (eq. 23). Panel a shows p{R)B? as a function of radius for 
a non-spinning BH. Panel h shows a similar plot for a* = 0.95. Panels c and d show disk 
thickness if as a function of radius for a=K = and 0.95, for disks with asymptotic values 
of H/R = 0.01 and 0.1. In panels a and b, the normalizations used for H/R = 0.01 and 
H/R = 0.1 are different. 
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Fig. 4. — Vertically integrated stress 2HaP{xlO'^) for a non-spinning BH. In both panels 
the standard disk model is plotted as a thick solid line. For the line types/colors defining 
the various models and the location of the ISCO, see Figure 1. (a) For H/R = 0.01, all four 
models are seen as indistinguishable from the standard model. (6) For the thicker disk, the 
models can be cleanly distinguished inside R ~ 15i?g. All numerical values correspond to 
G = c = M = l, M = l. 
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Fig. 5. — Similar to Figure 4, but for a^,=0.95. As in Figure 2, there is no variable-a model, 
and the ISCO is located a.t R = 1.937i?„. 
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Fig. 6. — Rate of energy dissipation RdL/dR as a function of radius R for models of a 
non-spinning BH (a and h) and a spinning BH (c and d). For the line types/colors defining 
the various models and the location of the ISCO, see Figures 1 and 2. The models shown for 
the thinner disks coincide with the standard disk model, whereas the thicker disks deviate 
somewhat from the standard model. 
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Fig. 7. — Spectra corresponding to the numerical disk models described in this paper for 
models of a non-spinning BH (a and b) and spinning BHs with thicker disks for a* = 0.8 (c) 
and a* = 0.95 (d). For the line types/colors defining the various models, see Figures 1 and 
2. Again, the different models are essentially indistinguishable in the case of the thin disk 
(a). 
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Fig. 8. — (a) Error in the spin estimate as a function of a for a BH with a true spin parameter 
of a* = 0. The quantity Aa* is equal to the value of a* obtained from fitting the model 
spectrum minus the true a*. The different curves correspond to different relative thicknesses 
HjR of the disk. (6) Error in the spin estimate for the variable-a profile as a function of 
disk thickness. In panel a, note the variable offset from zero error that occurs near ol = 
0.01, which is not visible in the previous figures for which the the different a models nearly 
coincide with the standard disk model. 
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Fig. 9. — (a) Error in the spin estimate Aa* for a non-spinning BH, calculated using spectral 
fitting and eq. (27) (solid lines) and from eq. (29) (dotted fines). The agreement is very good, 
showing that the simpler approach via eq. (29) is quite accurate, (b) Aa* as a function of a 
for BHs with different values of a*, calculated using equation (29). 
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Fig. 10. — The thick red sohd hne shows the vertically integrated stress profile as predicted by 
the standard disk model for a non-spinning BH (PW80 potential). The other lines show the 
stress profiles of three hydrodynamic disk models with advection parameter values / = 0.25 
(short-dashed line), / = 0.5 (dotted line) and / = 1 (long-dashed line). All three models use 
the variable a prescription (eq. 22), and their stress profiles have been scaled to match the 
standard model at i? = 15Rg. A logarithmic scale has been used to facilitate comparison 
with the MHD simulation result shown in Fig. 10 of HK02. 



